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Using the Jordan- Wigner fermionization in two dimensions we obtain the zz wave 
vector- and frequency-dependent structure factor for the spin-i isotropic XY model on 
a spatially anisotropic square lattice. We use the obtained results to discuss a role of the 
interchain interaction for the dynamic properties of quasi-one-dimensional systems. 
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Dynamic properties of the two-dimensional quantum spin models have attracted much attention for 
the last years mainly due to the discovery of high -temperature superconductivity of the layered copper 
oxides^. At present, many compounds are known to be good realizations of the square-lattice quantum 
spin model and the experimental investigations of their dynamic properties^ require corresponding the- 
oretical studies. In most cases the isotropic Heisenberg Hamiltonian is used to appropriately model the 
spin degrees of freedom. However, some generic features of the square-lattice quantum spin systems can 
be illustrated within the context of the spin-i isotropic XY model, which is far more amenable to a 



theoretical analysis. Moreover, a study of the latter model can be viewed as a first step for examining 
the more complicated case of the square-lattice spin-^ Heisenberg model. 

The aim of this study is to examine the dynamic properties of the spin-i isotropic XY model on 
a square lattice. Many theoretical and numerical investigations concerning the thermodynamic prop- 
erties of this model have been performed until now^. Less progress has been made in the studies of 
its dynamic properties. To examine the dynamic properties of the model we use the two-dimensional 
Jordan-Wigner fermionization^"'' (for a brief review see Ref. 8) thus extending the approach which was 
previously applied mainly to the study of thermodynamics (see also several recent papers^ on calcula- 
tion of the magnetization curves of some two-dimensional spin systems). More specifically, we calculate 
the two-site time-dependent correlation functions of z spin components and the corresponding dynamic 
structure factor ^'^^(k, w) for the spin-^ isotropic XY model on a spatially anisotropic square lattice. To 
trace a one-dimensional to two-dimensional crossover, we compare the results obtained with the corre- 
sponding ones for the spin-^ isotropic XY chain thus demonstrating the role of the interchain interaction 
for the dynamic properties of quasi-one-dimensional spin-^ isotropic XY systems. The recent study of 
the static properties of the square-lattice spin-| isotropic XY model^° suggests that the approach based 
on the two-dimensional Jordan-Wigner fermionization works well as long as the interchain interaction is 
small and it gives qualitatively correct results for larger values of the interchain interaction. 

We consider a spin model consisting of iV ^ oo spins ^ on a spatially anisotropic square lattice governed 
by the isotropic XY Hamiltonian 

oo oo 
i=0 j=0 

Here J > and J± > are the exchange interactions between the neighbouring sites in a row and 
in a column, respectively. We are interested in the two-site time-dependent spin correlation functions 
{^n,m{i)^n+p,m+q) (^^6 angular brackets denote the canonical thermodynamic average with the Hamilto- 
nian (1)) and the zz dynamic structure factor 

oo oo „QQ 

5,,(k,a;) = ^^ei('=^^+'='") / dte^-* («,„(*)<+,,„+,) - ■ (2) 

To calculate these quantities we use the two-dimensional Jordan-Wigner transformation^ reformulating 
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the problem in fermionic language. First we introduce the annihilation and creation operators of (spinless) 
fermion dij = e-^°''-isr., aij = T,%o{jii)T,'^=o{^j) BhrJ,9dtg^f,9' Bi,3;f,9 = Imln (/ - i + i (5 - j)). 
We adopt the mean-field treatment of the phase factors which appear in the initial Hamiltonian (1) after 
fermionization (for details see Refs. 5, 8). This is the only approximation made. Then we perform the 
Fourier transformation and the Bogolyubov transformation to arrive at the Hamiltonian of noninteracting 
fermions. Next, we rewrite the correlation function {Sn,mi^)^n+p,m+q) fermionic language and use the 
Wick-Bloch-de Dominicis theorem to calculate it. Finally, we insert the derived result into (2) and get 
the zz dynamic structure factor 

dkiy r dfci^ 
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cos^ 7ki ^^^^ _ ^^^^^-^ _ A^^_^^) 

+ (1 - rikj nki+ki^ (w - Aki + Aki+k)) 
+ sin2 'T^ki+k^- 7ki ^nk,nk,+i,6 {u + Au, + Ak^+k) 

+ (1 - nkj (1 - riki+k) S{u- Aki - Ak^+k))) • (3) 



Here Ak = ^ sin^ + Jj_ cos^ ky > is the fermion energy, cos ^ = \J\ + ~ 2Ak^" ' ^^"^ "2^ ~ 
sgn(Jsinfca;) \J\~ ~ 2Kk^ ' ''^^ ~ e^'^k+i Fermi factor. In Figs. 1 - 3 we present (k,a;) 

for the square-lattice spin-^ isotropic XY model as it follows from Eq. (3) (to integrate over k\x we 
use the relation 5{if{x)) = Y^j^ix — aj) \(p'{aj)\~^, aj are the roots of equation (p{x) = 0, and then 
we do the integral over kiy numerically). Let us discuss the dynamic properties of the considered spin 
model (1) concentrating in particular on the appearance of one-dimensional behaviour at low and high 
temperatures as ^ decreases. 

It is useful to start the analysis of dynamic properties by reminding what is long known for the one- 
dimensional case. From the exact calculation for spin-| isotropic XY chain^^ we know that Szzik,uj) 
at zero temperature /3 = 00 is governed by a two-fermion excitation continuum with lower and upper 
boundaries J| sinfc| and 2J| sin ||, respectively. The spectral weight is finite at the lower boundary. The 
lower boundary touches the line a; = at fc = 0, 27r and k = w (gapless modes). The spectral weight 
increases towards the upper boundary where Szz{k,uj) diverges. With the increasing temperature, the 



lower boundary gets smeared out and finally disappears, whereas the upper boundary remains unchanged. 
The described picture agrees with what we see from the greyscale plots of Szz{kx,Q, for ^ = 0.1 (Figs, 
la and 2a), i.e., for a system of almost noninteracting chains for momentum transfer along the chain 
direction. (The gap seen in Fig. 2a arises due to a finite interchain interaction: as J_l — > it shrinks 
and disappears.) For momentum transfer perpendicularly to the chain direction, the dynamic structure 
factor for a system of almost noninteracting chains disappears (Figs, lb and 2b). 

Let us pass to the two-dimensional case (the rest of the plots in Figs. 1, 2). Similar to the one- 
dimensional case SzzO^,i^) within the framework of our approach is conditioned by two-fermion excita- 
tions and hence it exhibits nonzero values only within the restricted region having a sharp high-frequency 
cutoff (evidently, S^z (k, w) (3) equals to zero if the frequency u) exceeds 2 \/ J'^ + J^) . Contrary to the one- 
dimensional case, for which Szzik^to) is almost structureless within the excitation continuum (apart from 
the upper boundary singularities), S'2z(k, w) exhibits several washed-out excitation branches (modes). 
Consider at first zero temperature. As it follows from Eq. (3), the (5-function 6 (uj — Ak^ — Ak^+k) admits 
the spectral weight at k, u) due to the two fermions with ki = ^— ^ — and ki + k = f + 
with the energy of the pair 

= 2^J2sin2^ + J^sin^^. (4) 

This excitation is known as the spin wave^. Note, that since the derivative of the argument of J-function 
with respect to k\x tends to zero (except for kx = 0, 2tt if ky = 0, 2tt and vise versa) the peaks of 
Szziy,!^) along the line uj\P (4) are expected. They are nicely visible in Figs. 1, 2. 

It is easy to note that the branch which almost everywhere in the regions {kx, 0, w) and (0, ky, to) forms 
the upper boundary (see Figs. 1, 2) is given by 

4!!o = 2/^W^7^, 4X=2^^^ + ^isin^^. (5) 

('2) ('2) 

^^^(k, w) arises along the lines and uJqI^ (5) when it is conditioned by the two fermions with 

(— O), (^, O) and § — , § + , respectively. There are other high-frequency excitations 
in Figs. 1, 2, conditioned by two fermions with ki = (0, 0) and ki + k {ky = 0) or with ki = (^, ^) and 
ki + k {kx = 0), i.e., 
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It is interesting to note, that the derivative of the argument of (5 (w — Ak^ — Ak^+k) with respect to kix 
tends to zero at fc^ = f , ^ {ky = 0) or at ky = ^, ^ {k^ = 0) that manifests itself as the increase of 
spectral weight in these regions clearly seen in some plots in Figs. 1,2. 

In the low-temperature limit, because of the Fermi factors, the ^-functions 6{ui — Aki + Ak^+k) and 
5{u) + Aki — Aki+k) in Eq. (3) may contribute only if the energy of one of the two fermions equals to 
zero, i.e., if either Aki = or Ak^+k = 0. As a result, we easily derive the lower boundary for the region 
of nonzero values of S^zi^,/^) 



At high temperatures, the Fermi factors do not impose the mentioned restriction. Putting kix = — 
kiy = ^ — kx ioT ky = and kix = ky, k\y = | — A;^ for fc^ = we get 



The excitation branch wj^ ' (8) contains most of the spectral weight at high temperatures (see Fig. 2). 

Obviously in our analysis of S^zi}^, w) we focus mainly on the consequences presumed by the ^-functions 
involved in (3) omitting the precise examination of a role of the rest of the factors (e.g., of the Fermi 
factors which controls a redistribution of the spectral weight as temperature increases). On the other 
hand, the results presented in Figs. 1, 2 (and Fig. 3) reflect all those underlying features of Eq. (3). 

In Fig. 3, we show the constant wave vector scans of zz dynamic structure factor, i.e., the frequency 
dependence of Szzi}^,^^) for certain wave vectors (the corresponding sections in Figs. 1, 2 can be easily 
located). All the discussed properties of ^^^(k, w) clearly manifest themselves in the depicted frequency 
profiles which may be almost symmetric or asymmetric, resemble (5-peaks, result from two coalesced 
peaks or exhibit tails which gradually disappear or are abruptly cut off. We indicate by numbers j some 
well pronounced peaks and cusps corresponding to the modes uj'^^ (4) - (8). 

To summarize, we have presented the first results for zz dynamic structure factor of the spin-^ isotropic 
XY model on a spatially anisotropic square lattice derived on the basis of the two-dimensional Jordan- 
Wigner fermionization. Our study has indicated excitations which govern the zz dynamic structure factor. 
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The established modes may manifest themselves as peaks or cusps in frequency profiles of (k,a;) and 
may be used for determining the Hamiltonian parameters. (For example, the gap seen in Fig. 2a and 
Fig. 3e is controlled by i<j'j^\, '^'k \^ ^^'^ '^k'\^ and, in particular, wlf g = ^y + J'^ + J±, wlf g = J, 
and wlfg = ^ J"^ -\- J\ — Jj_.) Thus, the theoretical result observed in our work should prove valuable 
in understanding the experimentally observable dynamic properties of two-dimensional spin-^ isotropic 
XY materials. However, there are, to our knowledge, no experimental results yet available (therefore, no 
compound specific parameters were considered above) that enable a direct comparison between theory 
and experiment to be made. 

It should be stressed that our results for the considered two-dimensional spin model contains the exact 
ones in the one-dimensional limit (after putting J_l = or J = 0). On the other hand, our approach is 
approximate because of the mean-field description of the phase factors which arise after fermionization. 
Such a treatment neglects the complicated interaction between fermions. It would be desirable to estimate 
the effects of this simplification going beyond the mean-field scheme for the phase factors. A possible 
test of the results for spin correlation functions is a comparison with the corresponding data of exact 
diagonalization and other numerical approaches. (Note, however, that the results reported in Refs. 3, 
10 refer to the in-plane spin correlation functions but not to the zz spin correlation functions which 
are different for the isotropic XY model in contrast to the case of the isotropic Heisenberg model.) In 
the case of the isotropic Heisenberg model the interaction between fermions is present even within the 
adopted approximate procedure due to the (Ising) interaction of z spin components. The quartic terms 
in the fermion Hamiltonian are treated after making further approximation^. One of the ways to proceed 
(the uniform flux solution) yields the same results for both the XY and Heisenberg models and from 
such a viewpoint our results for the dynamic properties may refer to the square-lattice spin-^ isotropic 
Heisenberg model as well. However, as one can conclude from the studies reported in Refs. 5, 6, a more 
sophisticated treatment of the fermion interaction inherent in the latter spin model (the in-phase flux 
solution, the Neel flux solution, the in-phase Neel flux solution) would be more appropriate. It seems 
obvious to extend the present investigation to the study of the square-lattice spin-i isotropic Heisenberg 
model which is used for interpreting the data for corresponding compounds^. It is interesting to note 



that the recent neutron scattering studies on CS2CUCI4 (the minimal model of this two-dimensional 
compound is the spin-5 antiferromagnetic Heisenberg model on a triangular lattice) definitely show a 
highly dispersive continuum of excited states. The conventional spin-wave theory clearly fails to account 
for the observed line shapes which do not show single particle poles but rather extended continua. We 
can also see some similarities at a qualitative level of the frequency profiles plotted in Fig. 3 of the present 
paper and in Fig. 3 of Ref. 12. However, we should bear in mind that CS2CUCI4 is a triangular but not 
a square-lattice two-dimensional system and that the minimal Hamiltonian determining the magnetic 
order is the isotropic Heisenberg model but not the isotropic XY model. New theoretical work is needed 
to elaborate the two-dimensional Jordan-Wigner fermionization results for this compound. Another 
possible application refers to dynamics of spin-| ladders. Recently the dynamic properties of the two- 
leg spin-^ ladder have been examined using the one-dimensional Jordan-Wigner fermionization. It seems 
interesting to compare these results with the corresponding ones which arise within the two-dimensional 
Jordan-Wigner fermionization approach. 

The authors thank S. Haas for a critical reading of the manuscript and helpful comments. O. D. is 
grateful to the Max-Planck-Institut fiir Physik komplexer Systeme (Dresden) for the hospitality during 
the International Workshop and Seminar on Modern Aspects of Quantum Impurity Physics (2003) when 
the paper was completed. 
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FIGURE CAPTURES 

FIGURE 1. The zz dynamic structure factor (greyscale plots) for the square-lattice spin-i isotropic 
XY model as it follows from Eq. (3) at zero temperature. J = 1, Jj_ = 0.1 (a, b), Jj_ = 0.5 (c, d), 
Jj_ = 0.9 (e, f); fcy = (a, c, e), fc^; = (b, d, f). 

FIGURE 2. The same as in Fig. 2 at high temperature /3 = 0.1. 

FIGURE 3. Frequency dependence of the zz dynamic structure factor (3) for momentum transfer along 
the chain {k^ = f (a, e), fc^; = tt (b, f), fcy = 0) and for momentum transfer perpendicularly to the chain 
{kx = 0, ky = ^ (c), ky = TT (d)) as the interchain interaction changes (J = 1, J_l = 0.1 (dotted curves), 
Jj_ = 0.5 (dashed curves), Jj_ = 0.9 (solid curves)) at zero temperature (a, b, c, d) and high temperature 
p = 0.1 (e, f). 
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